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Abstract. Fast thermalization, a subsequent strong buildup of elliptic flow of 
f~| . QCD matter and jet-quenching as found at the Relativistic Heavy Ion Collider 

Qh' (RHIC) are understood as the consequence of perturbative QCD (pQCD) in- 

teractions within the 3-f 1 dimensional parton cascade BAMPS. The main con- 
' tributions stem from pQCD inspired bremsstrahlung. Comparisons to Au+Au 

. data of the flow parameter V2 as a function of participation number as well as 

the gluonic contribution to the nuclear modification factor Raa for most cen- 
tral collisions are given. Also the shear viscosity to entropy ratio is dynamically 
J> ' extracted, which lies in the range of 0.08 and 0.15, depending on the chosen 

, coupling constant. 
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The values of the elliptic flow parameter V2 measured by the experiments at the 
Relativistic Heavy Ion Collider (RHIC) [[I] are (nearly) as large as those obtained 
from calculations employing ideal hydrodynamics. This finding suggests that a fast 
, local equilibration of quarks and gluons occurs at a very short time scale < 1 fm/c. 

^ ' and that the locally thermalized state of matter created, the quark gluon plasma 

(QGP), behaves as a nearly perfect fluid exhibiting strong "explosive" collective 
motion. Quarks and gluons should be rather strongly coupled, pointing towards a 
small viscosity to entropy coefficient for the QGP. Besides the strong collective flow, 
jet-quenching has been observed at RHIC [21 as the second striking new discovery. 
So far, both phenomena could not be related by a common understanding of the 
underlying microscopic processes: Jet-quenching is widely thought to be described 
by investigating the potential energy loss on the partonic level in terms of elastic 
pQCD collisions and pQCD radiative processes. On the other hand, the strong 
elliptic flow implies early thermalization and liquid like bulk properties which can 
not be understood by binary pQCD collisions being considerably too weak. This 
has raised the speculation about nonperturbative interactions as well as about super 
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symmetric representations of Yang-Mills theories using the AdS/CFT conjecture. 

In contrast, in this talk, it is demonstrated that perturbative QCD can still 
explain a fast thermalization of the initially nonthermal giuon system [ 131 131 [5] , 
the large collective effects of QGP created at RHIC [ [6] , the smallness of the shear 
viscosity to entropy ratio [ [7l [6j and also the jet-quenching of high momentum 
partons [ '8] in a consistent manner by using a relativistic pQCD based on-shell 
parton cascade Boltzmann approach of multiparton scatterings (BAMPS) [[3l [4]. 

BAMPS is a parton cascade, which solves the Boltzmann transport equation 
and can be applied to study, on a semi-classical level, the dynamics of gluon matter 
produced in heavy ion collisions at RHIC energies. The structure of BAMPS is based 
on the stochastic interpretation of the transition rate [ 3 , which ensures full detailed 
balance for multiple scatterings. BAMPS subdivides space into small cell units 
where the operations for transitions are performed. Gluon interactions included 
in BAMPS are elastic and screened Rutherford-like pQCD gg gg scatterings as 
well as pQCD inspired bremsstrahlung gg <-> ggg of Gunion-Bertsch type. The 
matrix elements are discussed in the literature [ [H HI [8] . The suppression of the 
bremsstrahlung due to the Landau-Pomeranchuk-Migdal (LPM) effect is taken into 
account within the Bethc-Heitler regime employing a step function in the infrared 
regime, allowing for independent gluon emissions. 

In the present simulations, the initial gluon distributions are taken in a Glauber 
geometry as an ensemble of minijets with transverse momenta greater than 1.4 GeV 
[Hj) produced via semihard nucleon-nuclcon collisions. The later interactions of the 
gluons are terminated when the local energy density drops below 1 GeV/fm'^. This 
value is assumed to be the critical value for the occurrence of hadronization, below 
which parton dynamics is not valid. Because hadronization and then hadronic cas- 
cade are not yet included in BAMPS, a gluon, which ceases to interact, propagates 
freely and can be regarded as a free pion employing a picture of parton-hadron du- 
ality. (Implementing a Cooper-Frye prescription for hadronization and employing 
UrQMD for the hadronic cascade are in progress.) The minijct initial conditions 
and the subsequent evolution using the present prescription of BAMPS for two sets 
of the coupling — 0.3 and 0.6 give nice agreements to the measured transverse 
energy per rapidity over all rapidities [ 6 . 

As a first example the fast thermalization of gluons is demonstrated in a local 
and central region which is taken as an expanding cylinder with a radius of 1.5 fm 
and within an interval of space time rapidity —0.2 < rj < 0.2. Figure [T] shows the 
varying transverse momentum spectrum at various early times obtained from the 
BAMPS calculations for central Au+Au collisions at = 200 AGeV. The coupling 
constant is taken as — 0.3. In contrast of employing only binary pQCD collisions, 
the spectrum reaches an exponential shape at 1 fm/c and becomes increasingly 
steeper at late times [ El H] . This is a clear indication for the achievement of local 
thermal equilibrium and the onset of (quasi-) hydrodynamical collective expansion 
with subsequent cooling by longitudinal work. 

As a note aside, we remark that, in addition, potential color instabilities [|9] may 
play a role in isotropization of particle degrees of freedom at the very initial stage 
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Fig. 1. (Color online) Transverse momentum spectrum in the central region at 
different times obtained from the BAMPS simulation of a central Au+Au collision. 



where the matter is super dense. However, more quantitative studies are needed to 
determine their significance on the true thermal equilibration as suggested for the 
expanding quark gluon matter at RHIC. 

The inelastic pQCD based bremsstrahlung and its back reaction are essential 
for the achievement of local thermal equilibrium at a short time scale. The fast ther- 
malization happens in a similar way if color glass condensate is chosen as the initial 
conditions. One of the important messages obtained there is that the hard gluons 
thermalize at the same time as the soft ones due to the ggg — > gg process, which is 
not included in the so called standard "Bottom Up" scenario of thermalization [[5]. 

Kinetic equilibration relates to momentum deflection. Large momentum deflec- 
tions due to large-angle scatterings will speed up kinetic equilibration enormously. 
Whereas the elastic pQCD scatterings favor small-angle collisions, the collision and 
emmission angles in bremsstrahlung processes are, for lower invariant, i.e. thermal 
energies, rather isotropically distributed due to the incorporation of the LPM cutoff 
[m 11]. Hence, although the elastic cross section is still considerably larger than the 
inelastic one, the given argument is the intuitive reason why the bremsstrahlung 
processes are acting more effectively in the equilibration of the gluon matter than 
the elastic interactions. Quantitatively it was shown in detail in [[3] that the con- 
tributions of the different processes to momentum isotropization are quantified by 
the so called transport rates 

r "^^p pj n _ I pL\ r jip_r' 
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where Ci[f] is the corresponding coUision term describing various interactions, i = 
gg — > gg,gg — ^ 999:999 ~^ 99i respectively. The sum of them gives exactly the 
inverse of the time scale of momentum isotropization, which also marks the time 
scale of overall thermalization. As it turns out, either by direct calculation or by 
the simulation, R^gg^ggg is a factor of 3 — 5 larger than R*Jg_,gg over a wide range in 
the coupling constant, which demonstrates the essential role of the bremsstrahlung 
for thermal equilibration [ 4 . For a gluon gas, which is initially far away from 
equilibrium, one can roughly estimate the time scale of thermalization r^q by taking 
the inverse of the sum of the transport collision rates close to thermal equilibrium. 
At temperature T = 400 MeV one has T^q^l/T, = 0.32 fm/c for Os = 0.3. We 
note that the above hinges on the assumption that the system is static. Expanding 
systems are more complicated because particles flow, which drives the system out 
of local equilibrium. Therefore, the momentum degradation of flowing particles 
toward isotropy is slower than the inverse of the total transport collision rate [3]. 
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Fig. 2. (Color online) Elliptic flow V2i\y\ < 1) from BAMPS using = 0.3, and 
0.6, compared with the PHOBOS [ IGJ and STAR [TT] data. 

Further, employing the Navier-Stokes approximation the shear viscosity 77 is 
directly related to the transport rate [[3, 

77 = —71 — 

5 E^*^ + |«9t(lnA) ' 

where A denotes the gluon fugacity. This expression allows to calculate the vis- 
cosity dynamically and locally in a full and microscopical simulation. On the 
other hand, close to thermal equilibrium the expression reduces to the more in- 
tuitive form r] — ^ and thus for the shear viscosity to entropy ratio ^ = 



■ STAR ¥^{4} 

• PHOBOS Track-based' 
^ PHOBOS Hit-based 




pQCD description of elliptic flow and jet-quenching 



5 



{pl3R}gg_^gg + ^ ^ ^ Q Q ' Within thc present description bremsstrahlung and 
its back reaction lower the shear viscosity to entropy density ratio significantly by 
a factor of 7, compared with the ratio when only elastic collisions are considered. 
For tts = 0.3 one finds rj/s = 0.13 [[71 [5]. To match the lower bomid of ij/s = l/Air 
from the AdS/CFT conjectm'e = 0.6 has to be employed. Even for that case the 
cross sections are in the order of 1 mb for a temperature of 400 MeV. Perturbative 
QCD interactions can drive the gluon matter to a strongly coupled system with an 
rj/s ratio as small as the lower bound from the AdS/CFT conjecture. 

The elliptic flow V2 can be directly calculated from microscopic simulations and 
compared to experimental data, assuming parton-hadron duality, as = 0.3 and 0.6 
are used for comparisons [ 6]. Inspecting Figure [21 except for the central centrality 
region the results with — 0.6 agree perfectly with the experimental data, whereas 
the results with ~ 0.3 are roughly 20% smaller. One can summarize that suffi- 
cient elliptic flow is being built up, yet the details of the freeze-out, hadronization 
and possible hadronic phase contributions might also give some minor dependence 
on the overall strength in elliptic flow, which still have to be addressed. In any case, 
for the QGP phase gluon bremsstrahlung dominates and yields rapid thermaliza- 
tion, and, therefore, early pressure buildup, and a small shear viscosity like a nearly 
ideal fluid. 
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Fig. 3. (Color online) Shear viscosity to entropy density ratio rj/s at the central 
region during the entire expansion, rj/s values are extracted from the simulations 
at impact parameter 6 = 4, 6.3, and 8.6 fm. The upper band shows the results with 
as = 0.3 and the lower band the results with as = 0.6. 

Within these simulations the ratio of the shear viscosity to the entropy density, 
rj/s, can be locally extracted by means of the above formula. This is shown in Figure 
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[21 As expected from the explicit expressions of the matrix elements, the ratio does 
not depend strongly on the gluon density or temperature, since interaction rates and 
transport collision rates scale with the temperature. Hence, rj/ s depends practically 
only on Ug- For = 0.6, at which the V2 values match the experimental data, one 
has rj/ s ~ 0.08. However, 77/s may be higher, since inclusion of hadronization 
and subsequent hadronic cascade can yield minor contributions to the final elliptic 
flow values. Furthermore, different picture of initial conditions (eg. color glass 
condensate) will also lead to different initially spatial eccentricity, and, hence, will 
affect the final value of V2 and the result on rj/s. These investigations are underway 
and will provide more constraints on extracting rj/s. 
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Fig. 4. (Color online) Gluonic Raa at midrapidity {y e [—0.5, 0.5]) as extracted from 
simulations for central Au+Au collisions at 200 AGeV. For direct comparison the 
result from Wicks et al. [[ij for the gluonic contribution to Raa and experimental 
results from PHENIX [[12] for 7r° and STAR [[T3] for charged hadrons are shown. 

The quenching of gluonic jets can now self-consistently be addressed within 
the same full simulation of BAMPS [ [8] including elastic and radiative collisions. 
Due to the steeply falling spectrum of initial mini-jets a suitable weighting and 
reconstruction scheme has to be employed, where a huge number of initial spectra 
are sampled which are then selected for such events for further simulation that 
contain high-pT^ partons. The nuclear modification factor Raa of the gluons is 
obtained directly by taking the ratio of the final px spectra to the initial mini-jet 
spectra. Fig. [ijshows the result, exhibiting a clear suppression of high-py gluon jets 
at a roughly constant level of R^a"^^ « O.f , potentially slightly rising towards high 
Pt- The coupling constant is taken as — 0.3. The dominant contribution for the 
energy loss are the bremsstrahlung contributions, although the picture changes for 
higher jet energies [[8]. The suppression is approximately a factor of two stronger 
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than the experimental pion data. This, however, was to be expected since at present 
the simulation does not include quarks, which are expected to lose less energy by a 
factor of 4/9. Indeed, comparing with state of the art results from Wicks et al. [ 
[14] for the gluonic contribution to Raa (seen as the line in Fig. |4]), which in their 
approach together with the quark contribution reproduces the experimental data, 
one finds a perfect agreement. 

In summary, the pQCD based parton cascade BAMPS is used to calculate the 
time scale of thermalization, the elliptic flow V2, extraction rj/s, and the quenching 
of gluonic jets from simulations of Au+Au collisions at RHIC energy ^/s = 200 
AGeV within one common setup. This is a committed and large scale undertaking. 
BAMPS includes elastic gg — > gg and inelastic bremsstrahlung and its back reaction 
99 999- The present approach thus constitutes a flrst realistic 3+1 dim. micro- 
scopic transport simulation that can describe the two phenomena elliptic flow and 
jet-quenching by the same underlying processes. Further analyses on jet quenching, 
on particle correlations, on quark degrees of freedom including hadronization, on 
initial conditions and on the exploration and use of dissipative hydrodynamics are 
underway to establish a more global picture of heavy ion collisions. 
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